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Abstract. Central supermassive black holes (SMBHs) are a ubiquitous feature of 
locally-observed galaxies, and ample evidence suggests that the growth of SMBHs and 
their host galaxies is closely linked. However, in the event of a merger, gravitational- 
wave (GW) recoil may displace a SMBH from its galactic center, or eject it entirely. 
To explore the consequences of this phenomenon, we use hydrodynamic simulations of 
gaseous galaxy mergers that include a range of BH recoil velocities. We have generated 
a suite of over 200 simulations with more than 60 merger models, enabling us to iden- 
tify systematic trends in the behavior of recoiling BHs - specifically (i) their dynamics, 

(ii) their observable signatures, and (iii) their effects on BH/galaxy co-evolution, (i) Re- 
coiling BH trajectories depend heavily on the gas content of the host galaxy; maximal 
BH displacements from the center may vary by up to an order of magnitude between 
gas-rich and gas-poor mergers. In some cases, recoil trajectories also depend on the 
timing of the BH merger relative to the formation of the galaxy merger remnant, (ii) 
Recoiling BHs may be observable as offset active galactic nuclei (AGN) via either kine- 
matic offsets (v > 800 km s"') or spatial offsets (R > I kpc) for lifetimes of about 1 - 
100 Myr In addition, recoil events affect the total AGN lifetime. GW recoil generally 
reduces the lifetimes of bright AGN, but may extend lower-luminosity AGN lifetimes. 

(iii) Rapidly-recoiling BHs may be up to about 5 times less massive than their station- 
ary counterparts. These mass deficits lower the normalization of the M - cr relation 
and contribute to both intrinsic and overall scatter Furthermore, recoil events displace 
AGN feedback from the galactic center, which enhances central star formation rates. 
This results in longer starburst phases and higher central stellar densities in merger 
remnants. 



1. Introduction 



Recent results from numerical relativity simulations indicate that "central" supermas- 
sive black holes (SMBHs) may in fact spend substantial time in motion and offset from 
the galactic center, owing to the large gravitational-wave (GW) recoil kicks that may 
accompany BH mergers (e.g., |Campanelli et al.|2007[ [Baker et aL]|2008 Lousto et al. 



20101. RecoiUng BHs m ay be obs ervable as offset AGN ([Madau & Quataert||2004 
Loeb|2007HBlecha & L oeb 2008; Bl echa et al.|20TT]|Guedes et al.|20imSijacki et al. 



201 1 1. Several such candidates have been discovered, though none have yet been con- 
firmed ([Koinossa et al.|2008t|Shields|2009t|Civano.2010 ; Batcheldor et al.|20T0t|Jonker| 



et al.||2010 i. There may also be indirect consequences of GW recoil for BH 
co-evolution ( |Volonteri|2007| [Blecha et al.|2011[ [Sijacki et al.|2011 1. 
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Figure 1. Recoil trajectories in gas-rich (left panel, 60% gas initially) and gas- 
poor (right panel, 10% gas) mergers. In each case, the recoiling BH separation from 
the galactic center is plotted versus time for varying Vk/Vesc within a single galaxy 
model, v'k/vesc = 0.6, 0.8, & 0.9 are shown, as indicated on the plots. 



We conduct a parameter study of recoiling BHs, using hydrodynamic simulations 
of galaxy mergers. Our simulations are described in § |2] We discuss recoiling BH 
dynamics in § |3.1| In § |3.2[ we examine the offset lifetimes of recoiling AGN, and in 



§ |3.3[ we examine the effects on host galaxies. Our conclusions are summarized in § ffl 



2. Galaxy Merger Simulations 



We simulate galaxy mergers using the smoothed particle hydrodynamics (SPH) code 
GADGET-3 ( |Springel|2005[ [Springe! et al.|2005) . We add an arbitrary kick velocity Vk to 
the remnant BH at the time of BH merger (fmrg). which we generally scale to the central 
escape speed, Vesc(fmrg)- The BHs are allowed to accrete both from ambient gas via the 
Bondi-Hoyle formula and from an ejected disk of gas with a time-dependent accretion 
rate. These prescriptions are described in detail in Blecha et al. ( 2011| l. 

We have used 62 different galaxy merger models in which we vary the galaxy mass 
ratio (q), the total galaxy mass, the gas fraction (/gas), and the orbital configuration. For 
each model, we simulate both a merger with no recoil kick and a merger with Vk/vesc - 
0.9. For a subset of these models, we also simulate intermediate values of v^. We refer 
the reader to [Blecha et al., ( ,2011| ) for the full details of our initial conditions. 



3. Results 

3.1. Recoiling BH Trajectories 

Fig.[T]demonstrates that the gas content of galaxies greatly influences recoiling BH tra- 
jectories. In the left panel (right panel), trajectories with = 0.6, 0.8, &0.9 Vgsc are 
shown for a merger in which the progenitor galaxies had 60% (10%) gas initially. In 
both cases, BHs kicked with Vk = 0.9 Vesc are still on large orbits by the end of the sim- 
ulation. However, at more moderate kick speeds, the BH has a much shorter wandering 
time in the gas-rich remnant than in the gas-poor remnant. For Vk = 0.6 Vesc, the BH in 
the gas-rich remnant is confined to the central kpc. 



Recoiling Black Holes in Merging Galaxies 
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Figure 2. Left panel: Kinematically-ofFset AGN lifetimes for varying Vk/vgsc with 
^min = 500 km s \ for the 3 AGN definitions given in the text. Middle panel: same 
for Vmin - 800 km s Right panel: spatially-offset AGN lifetimes fori^niin - 1 kpc. 



3.2. AGN Lifetimes 

We calculate the observable lifetimes of recoiling BHs as spatially- or kinematically- 
offset AGN (Jagn)- We use 3 different definitions of an AGN: L^oi ^ (10%^Edd> 
3% LEdd> 3.3 X 10^ L0). Lboi & L^dd are the bolometric and Eddington luminosities. 
For kinematically (spatially) offset AGN, we require a BH velocity offset vbh > ^min 
(^BH > ^min) from the Stellar center of mass. Fig. |2] shows the resulting offset AGN 
lifetimes. At moderate kick speeds, a kinematically-offset AGN may occur when the 
BH passes through the gas-rich galactic center. At higher kick speeds with longer or- 
bits, the BH appears as a kinematically offset AGN only immediately after the kick. 
Spatially-offset AGN may have longer lifetimes, but generally at lower luminosities. 

3.3. Co-evolution of BHs and Galaxies 

Recoiling BHs have smaller final masses than stationary BHs, in some cases by up to a 
factor of ~ 5. In the left panel of Fig.|3} we compare the Mbh - cr* relation resulting 
from our sets of simulations with no recoil and with large recoil (vk/vesc = 0.9). These 



data do not reproduce the observed Mbh - c* relation (e.g., Tremaine 2002 1, as we have 



varied parameters systematically, but the relative differences between the simulated 
correlations are informative. GW recoil lowers the offset of the Mbh - o"* relation and 
increases intrinsic scatter. In a cosmological framework, scatter would increase even 
further as some recoiling BHs were replaced by new BHs via subsequent mergers. 

GW recoil can also affect the central structure of galaxy merger remnants. The 
right panel of Fig.[3]compares the star formation rates (SFRs) for simulations of a gas- 
rich merger with v^/vesc - and 0.9, and with no BHs. The simulation with large recoil 
has a higher post-merger SFR than the no-recoil simulation, similar to the case with no 
BHs. This is a direct result of the displacement of AGN feedback via GW recoil. 



4. Conclusions 

We find that recoiling BHs may appear as offset AGN for ~ 1-100 Myr, for a wide 
range in kick speeds. Their trajectories are greatly influenced by the environment of a 
gaseous merger remnant. Moreover, recoiling BHs themselves may influence the cen- 
tral regions of the merger remnant and contribute to scatter in the BH - bulge relations. 



4 



Blecba, Cox, Loeb, and Hemquist 




-0.6 -0.4 -0.2 0.0 0.2 0.4 1 2 3 4 

log(<7 / 200 km s') ^ime [Gyr] 



Figure 3. Left: Mbh - o", relation for simulations with no recoil kicks (open 
circles) and with Vk/^esc = 0.9 (filled triangles). Mbh is the final BH mass and cr, is 
the line-of-sight-averaged stellar velocity dispersion; error bars show the range of o-, . 
The solid & dashed lines are least-squares fits to the data. Right: Total SFR vs. time 
for a gas-rich merger The simulations shown have Vk = (dashed line), Vk/vesc = 0.9 
(solid line), and no BHs (dotted line). The vertical line marks the time of BH merger 



References 

Baker, J. G., Boggs, W. D., Centrella, J., Kell y, B. J. McWilHams, S. T., Miller, M. C, & van 

Meter, J. R. 2008, ApJ, 682, L29. (5802.Q416| 
Batche ldor, P., Robinson, A., Axon, D. J., Perlman, E. S., & Merritt, D. 2010, ApJ, 717, L6. 

I1QQ5.2173 

Blecha, L., Cox, T. J., Loeb, A., & Hemquist, L. 2011, MNRAS, 412, 2154. |1089 ■ 494Q| 
Blecha, L., & Loeb, A. 2008, MNRAS, 390, 1311. 08(95 . 1420 



CampanelH, M., Lousto, C, Zlochower, Y., & Merritt, D. 2007, ApJ, 659, L5. |arXiv: gr-qc/ 
0701164 

Civano, F. effl/.. 2010, ApJ, 717, 209. 1003.00201 

Guedes, J., Madau, R, Mayer, L., & Callegari, S. 201 1, ApJ, 729, 125. 1008 . 2032 

Jonker, R G., ToiTes, M. A. R, Fabian, A. C, Heida, M., Miniutti, G., & Pooley, D. 2010, 

MNRAS,407, 645. 1004.5379 

Komossa, S., Zhou, H., & Lu, H. 2008, ApJ, 678, L81. arXiv 10804.4585 

Loeb, A. 2007, Physical Review Letters, 99, 041103. arXiv:astro-ph/0703722 

Lousto, C. O., Campanelli, M., Zlochower, Y., & Nakano, H. 2010, Classical and Quantum 

Gravity, 27, 1 14006. 0904 .3541 
Madau, P, & Quataert, E. 2004, ApJ, 606, L1 7.|arXiv:astro-ph/0403295| 

Shields, G. A. et al. 2009, ApJ, 707, 936. 0907 . 3470 

Sijacki, D., Springel, V., & Haehnelt, M. G. 2011, MNRAS, 414, 3656. 1 1008 ■ 31T3 
Springel, V. 2005, MNRAS, 364, 1105. arXiv: astro-ph/0505010 



Springel, V., Di Matteo, T., & Hernquist, L. 2005, MNRAS, 361, 776. |arXiv:astro-ph/ 
0411108 

Tremaine, S. et al. 2002, ApJ, 574, 740. arXiv:astro-ph/020 3468l 
Volonteri, M. 2007, ApJ, 663, L5. arXiv:astro-ph/0703180 



